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Abstract. We present new X-ray observations of the X- 

' ray pulsar IE 2259+586, obtained during March 1997, 

I with the Rossi X-Ray Timing Explorer (RXTE). We have 

■ measured the pulse frequency derivative v — (—1.08 ± 

\ 0.04) X 10~^^ Hz s~^ from pulse arrival times obtained 

' in a sequence of 5 observations spread over one month. 

' This v is consistent with the long term spin-down trend. 

' We also found that the observed X-ray luminosity is con- 

\ sistent with that measured at quiescent X-ray flux lev- 

' els by previous missions. Our observations imply that IE 

\ 2259-1-586 was spinning down steadily without exhibiting 

' any stochastic torque noise fluctuations during the month 

' covered by our observations. 
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1. Introduction 

The X-ray pulsar IE 2259-1-586 is located at the geo- 
metric center of the semi-circular shell of diffuse X-ray 
emission of the supernova remnant G109— 1.0 (Gregory & 
Fahlman 1980). Radio observations indicate that the age 
of G109— 1.0 is about 10^ yr and the distance to it is 3.6- 
4.7 kpc (Hughes et al. 1984). The pulse frequency of the 
pulsar is 0.1433 Hz (pulse period 6.8 sec) and the average 
spin down rate since its discovery is v = —1 x 10"^** Hz 
s~^. If IE 2259-1-586 is an isolated, spinning-down pulsar, 
this rate yields a rotational energy loss of 10'^"'^ erg 
which is much less then the intrinsic X-ray luminosity of 
~1035 erg s-i (Corbet 1995). Koyama et al., (1987, 1989) 
claimed that the spin-down rate and unabsorbed X-ray 
luminosity can be understood if IE 2259-1-586 is in a bi- 
nary system with a neutron star magnetic field of 5x10^^ 
Gauss and spinning close to the equilibrium period (Ghosh 
& Lamb 1979). However, efforts to find orbitally induced 
Doppler shifts in the X-ray pulsations which would reveal 
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the binary period have been unsuccessful. The upper lim- 
its of ax/csini from previous missions such as Einstein, 
EXOSAT, Ginga and ROSAT have been very small. For 
example, recent RXTE observations (Mereghetti et al., 
1998) have given the lowest upper limit ax/csini < 30 
msec at the 99% confidence level. According to these new 
upper limits, if IE 2259-1-586 is in a binary system, its 
companion star must be either a white dwarf, or a helium- 
burning star with M < O.SMo (Mereghetti et al. 1998). Al- 
ternatively, IE 2259-1-586 could be an isolated star which 
is accreting from a disk and is formed by remnants of the 
common envelope evolution of a high-mass X-ray binary 
(van Paradijs et al. 1995, Ghosh et al. 1997). 

Analyses of the ROSAT observations (Baykal & Swank 
1996) have shown that the source is not steadily spinning- 
down. A spin-up episode which is superposed on a long 
term spin-down trend was found. The deviations from the 
secular spin-down trend are also consistent with other ac- 
cretion powered X-ray pulsars (Baykal & Ogelman 1993), 
which are several orders of magnitude greater than those 
of radio pulsars. This result favors accretion onto a neu- 
tron star from a very low mass companion or from a resid- 
ual disk. In this research note, we present the pulse timing 
resuhs and X-ray spectrum of IE 2259-1-586 from RXTE 
observations. The main goal of these new RXTE observa- 
tions is to resolve the pulse frequency derivative (ly) during 
the observation time span and measure the source X-ray 
luminosity. These quantities enable us to further test the 
current accretion hypothesis for this source. 



2. Observation and Data Analysis 

IE 2259-1-586 was observed between 1997 February 25 and 
March 26 on 5 different days separated from each other 
by approximately one week. Each 25 ksec observation was 
taken in approximately one day. The results presented 
here are based on data collected with the Proportional 
Counter Array (PGA, Jahoda et al. 1996). The PC A in- 
strument consists of an array of 5 proportional counters 
operating in the 2-60 keV energy range, with a total effec- 
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tive area of approximately 7000 cm^ and a field of view, 
~ 1° FWHM. 

2.1. Torque and X-ray Luminosity 

Background light curves and the pulse height amplitudes 
are generated using background estimator models based 
upon the rate of very large events (VLE) , spacecraft acti- 
vation and cosmic X-ray emission with the standard PCA 
analysis tools and are subtracted from the source light 
curve obtained from the Good Xenon event data. After 
correcting the background subtracted light curves with 
respect to the barycenter of the solar system, data sets 
in each observation were folded on statistically indepen- 
dent trial periods (Leahy et al. 1983). A master pulse was 
constructed by folding the data on the period giving the 
maximum The master pulse with 55 phase bins was 
represented by its Fourier harmonics (Deeter & Boynton 
1985) and cross-correlated with the harmonic representa- 
tion of average pulse profiles from each observation. The 
pulse arrival times so obtained are represented in Fig. 1. 
The quadratic trend of pulse arrival times is a direct mea- 
sure of the change of the pulse frequency during the ob- 
servation (or intrinsic pulse frequency derivative). 



ho + 6u{t - to) + -u{t - tof 



(1) 



where dcf) is the pulse phase offset deduced from the pulse 
timing analysis, to is the mid-time of the observation, 
(j)o is the phase offset at to, 5y is the deviation from 
the mean pulse frequency (or additive correction to the 
pulse frequency), and v is the pulse frequency derivative 
of the source. Table 1 presents the timing solution of IE 
2259-1-586 from our RXTE observations. The pulse fre- 
quency derivative obtained from the quadratic trend of the 
pulse timing analysis is vrxte = —(1.08 ± 0.04) x 10"^'* 
Hz s~^ which is consistent with the average spin-down rate 
over a time span of 19 years, < v >= — (1.15±0.06) x 10^^^ 
Hz s^^. Fig. 2 presents the previous pulse frequency his- 
tory of IE 2259-1-586 together with the new RXTE obser- 
vation. 



Table 1. Timing Solution of IE 2259-1-586 for RXTE Obser- 
vations 




Fig. 1. Phase ofsets in pulse arrival times. 
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Fig. 2. Pulse frequency history of 1E2259-I-586, including the 
latest RXTE measurement. 



Epoch(MJD) 


50516.494628 


Pulse Frequency (Hz) 
Pulse Preq. Derivative (Hz s~^) 


0.1432886474 ± 1.8x10"'' 
-(1.08 ±0.04) X 10"" 




channel energy (keV) 



The background spectrum was calculated using the 
same background estimator models based upon the rate 
of very large events (VLE), spacecraft activation and cos- 
mic X-ray emission as used to calculate background light 



Fig. 3. (top) RXTE PCA X-ray spectrum of 1E2259-I-586 with 
the fit to the power law model, (bottom) The residuals to the 
fit. The reduced = 11- 
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Fig. 4. Results of a search for orbital Doppler delays for IE 
2259+586. The horizontal lino denotes the a;r/c sin i values 
corresponding to detection at the 99% confidence level. 

curves. The resultant background subtracted spectrum 
in the energy range 2-10 keV yields X-ray luminosity 
= (0.38±0.06) X 10^^(41^)^ erg s'^ which is consistent 
with that obtained from ASCA and SAX measurements 
in 2-10 keV (Corbet et al., 1995, Parmar et al., 1997). The 
best fitting power law X-ray spectrum with column den- 
sity Nh=(2.2± 0.8) X 10^^ cm-2 and photon index 4.78 is 
presented in Fig. 3. The background contamination from 
the supernova remnant is only a few percent of the total 
source flux (Rho & Petre 1997). ASCA and Beppo SAX 
observations have found a soft black-body component with 
a temperature of 0.44 keV (Corbet et al. 1995, Parmar et 
al. 1997). Our fits to the X-ray spectra do not resolve 
the soft black-body component since the RXTE/PCA de- 
tectors are insensitive to energies less than about 2 keV. 
There is no evidence for any deviation from a power law 
that might be attributable to a cyclotron feature. 

2.2. Search for Orbital Period 

In searching for orbital Doppler modulation, orbital pe- 
riods from ~ 500 sec to ~ 2.5 x 10^ sec were considered 
by oversampling the independent Fourier periods (Leahy 
1983) by a factor of 4. For each trial orbital period, 8 
pulse profiles are obtained at the different intervals of the 
orbital phase by correcting the light curve for the secular 
spin-down rate and folding the light curve at the the pulse 
period. The resulting pulse profiles were cross-correlated 
with the master pulse, which was obtained from all obser- 
vations. To search for any sinusoidal modulation in the 8 
pulse arrival times, Fourier amplitudes of the arrival times 
were computed and are presented as a function of the trial 
orbital period in Fig. 4. In the absence of an orbital signal, 
the expected distribution of Fourier amplitudes is the ^ 
distribution with 2 degrees of freedom (van der Klis 1989). 
As seen in Fig. 4, no significant signal detection is found 



at the 99% confidence level. The upper limit for any unde- 
tected signal is estimated as ax smi < 0.028 light-s, which 
is very close to the earlier upper limit given by Mereghetti 
et al. (1998). (We note that a slight improvement of the 
upper limit arises from the better statistics due to the 
longer RXTE observation with respect to previous obser- 
vations) The above upper limits on Doppler shifts confirms 
that if the IE 2259-1-586 is in a binary system, then the 
companion star must be either a white dwarf, or a helium- 
burning star (Mereghetti et al., 1998). 

2.3. Torque and X-ray Luminosity Changes in Accretion 
Powered Disk Fed Binaries and IE 2259+586 

In disk- fed systems, the accretion onto the neutron star 

is believed to be from a Keplerian disk (Ghosh & Lamb 
1979). The torque on the neutron star is given by 

itl = n{ws)M Ik, (2) 

where Ik = {GMroY^"^ is the specific angular momentum 
added by a Keplerian disk to the neutron star at the inner 
disk edge « 0.5rA, where = (2G'M)-i/V^/^M-2/7 
is the Alfven radius, fi is the neutron star magnetic mo- 
ment, n{ws) ~ 1.4(1 — Ws/'Wc)/{l — Ws) is a dimen- 
sionless function that measures the variation of the ac- 
cretion torque as estimated by the fastness parameter 

= n/nK{ro) = 27rP-iG-i/2^-5/V/^Ai"-3/7. Here 
Wc is the critical fastness parameter at which the accre- 
tion torque is expected to vanish {wc ^ 0.35 — 0.85 de- 
pending on the electrodynamics of the disk, Lamb 1989). 
In this model, the torque will cause a spin-up if the neu- 
tron star is rotating slowly {wg < Wc) in the same sense 
as the circulation in the disk, or spin-down, if it is rotating 
in the opposite sense (see Lamb 1991). Even if the neu- 
tron star is rotating in the same sense as the disk flow, 
the torque will spin-down the neutron star if it is rotat- 
ing too rapidly {wg » iWc)- In such a model one should 
sec positive correlation between angular acceleration (Q) 
and mass accretion rate (M) if the disk is rotating in the 
same sense as the neutron star. If the flow is from Roche 
Lobe overflow then the accreting material carries positive 
speciflc angular momentum I, therefore it is hard to imag- 
ine accretion flow reversals and hence the spin-up/down 
torques should be correlated with mass accretion rate M . 

In Ginga observations the source had flux levels a fac- 
tor of two higher than average (Iwasawa et al. 1992). 
This implied that the mass accretion rate on the source 
is indeed variable and fluctuations in the spin-down rate 
and possible spin-up trends should be expected (Baykal 
& Swank 1996). However, our RXTE observations found 
the source with an X-ray flux consistent with the secular 
spin-down. Therefore, it is quite natural to find the secu- 
lar spin-down rate from pulse arrival times, although this 
observation gives the spin-down rate of the source over 
the shortest observing interval since the discovery of the 
source. 
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Recent observations of accreting neutron stars have 
shown stochastic spin- up/down trends on time scales from 
days to a few years (Bildsten et al. 1997). In intervals 
between stochastic changes disk-fed sources show secular 
spin-up or down trends with lower values of noise strength 
(Baykal 1997). Some of the sources switch from spin- up to 
spin-down states without showing great changes in their 
mass accretion rates (Bildsten et al. 1997). These unusual 
behaviors led Baykal (1997) and Nelson et al. (1997) to 
the possibility of retrograde circulation of accretion disks. 
GX 4+1 shows correlation between the X-ray flux and the 
spin-down rate (Chakrabarty et al. 1997) which may sug- 
gest a retrograde accretion disk. In the application of the 
above observational results and current accretion theory 
to IE 2259-1-586, the spin-down rate should decrease or it 
should switch to spin-up at higher accretion rates (higher 
X-ray flux). On the other hand if the source has a very 
unusual counter-rotating accretion disk, in the high state 
the spin-down rate should increase with respect to the low 
state. Future X-ray observations at a high state could ad- 
dress this important question. 
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